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Antiproton-nucleus electromagnetic annihilation as a way to 
access the proton timelike form factors 
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Abstract Contrary to the reaction p p — > e + e~ with a high momentum incident antiproton on a free target 
proton at rest, in which the invariant mass M of the e + e~ pair is necessarily much larger than the pp mass 
2m, in the reaction pd — > e + e~n the value of M can take values near or below the pp mass. In the antiproton- 
deuteron electromagnetic annihilation, this allows to access the proton electromagnetic form factors in the 
time-like region of q 2 near the pp threshold. We estimate the cross section dap^ e + e - n / dM for an antiproton 
beam momentum of 1.5 GeV/c. We find that near the pp threshold this cross section is about 1 pb/MeV. 
The case of heavy nuclei target is also discussed. Elements of experimental feasibility are presented for the 
process pd — > e + e~n in the context of the PANDA project. 
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1 Introduction 

A lot of efforts are devoted to understanding the struc- 
ture of the nucleon, the building block of matter. The un- 
derlying theory is the field theory of strong interaction, 
QCD, in its non-perturbative regime. Electromagnetic pro- 
perties of the nucleon are fundamental pieces to this puz- 
zle. Among them, the electromagnetic form factors of the 
proton and the neutron are basic observables, which are 
the goal of extensive measurements. In the spacelike re- 
gion, i.e. for a virtual photon four-momentum squared 
q 2 < 0, these form factors give information about the spa- 
tial distribution of electric charge and magnetization in- 
side the nucleon. In the timelike region (q 2 > 0) they tell 
us about the dynamics of the nucleon-antinucleon (NN) 
interaction. 

A fully consistent description of the nucleon form fac- 
tors should include both domains, of spacelike and time- 
like q 2 , since these domains are related by crossing sym- 
metry. Such theoretical models are generally based on 
dispersion relations [ 1 , 2 , 3] or semi-phenomenological ap- 
proaches |U|5l. They predict a smooth behavior of the 
form factor in the measured regions, but a peaked be- 
havior in the timelike region below the NN threshold 
(0 < q 2 < Am 2 , where m is the nucleon mass), due to 



poles in the amplitude (see e.g. fig. [TJ taken from [6]). 
These poles are phenomenological inputs, built from me- 
son exchange, and their properties are fitted to the data 
in the measured regions. The corresponding irregulari- 
ties in form factors are related to the transition of pp to 
vector mesons which can decay in e + e~ pair via a virtual 
photon. 

The mesons with a mass near the pp mass can have 
a quasinuclear nature, i.e., they can be formed by bound 
states and resonances in the pp system. Such vector mesons 
were predicted in the papers HZUSl. Note that such mesons 
can be formed not only in the pp system but in NN in 
general and they can have not only vector quantum num- 
bers. A review on quasinuclear mesons in the NN system 
is given in [9]. 

The under-threshold region (0 < q 2 < 4m 2 ) is called 
unphysical because it cannot be accessed experimentally 
by an on-shell process. Some experiments have been per- 
formed in the vicinity of the NN threshold, either in pp — ► 
e + e~ at LEAR [10] or in the inverse channel e + e~ — > pp at 
Babar ITU , but they cannot go below this physical thresh- 
old. However, a nucleus provides nucleons with various 
momenta, in modulus and direction, and also various de- 
grees of off-shellness. Therefore it offers the possibility 
to produce an NN electromagnetic annihilation with an 
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Figure 1. Experimental data and predictions for the magnetic 
proton form factor in the domain —10 GeV 2 /c 2 < q 2 < 
10 GeV 2 /c 2 . The fi eure is taken from 161 . 

invariant mass squared q 2 = Sp^ smaller than 4m 2 . The 
main purpose of our paper is to explore this possibility, 
which may give access to the proton form factors in the 
underthreshold region, for an off-shell nucleon. The idea 
to use a nucleus for that purpose was explored in the 80's 
using deuterium (12). The reaction is then: 

pd — > e + e~n (1) 

(a crossed-channel of deuteron electrodisintegration). The 
aim of the present paper is to revive this study in view of 
the future antiproton facility FAIR at GSI. 

Other channels can give access to the off -shell nucleon 
form factors in the timelike region, including the under- 
threshold region; such processes have been studied the- 
oretically in ref. p"3) (jp — > pe + e~) and in refs. [14,15| 
(pp — > 7r°e + e~). 

The paper is organized as follows: a theoretical study 
is presented in sect. [2] and experimental aspects are pre- 
sented in sect. |3] Other aspects are mentioned in sect. H] 
and a conclusion is given in sect. [5] 

2 Theoretical study 

In elastic electron scattering from the nucleon e~N — > 
e~N the momentum transfer squared q 2 = (ft — ft:') 2 is al- 
ways negative. This allows to measure the nucleon form 
factors in the space-like domain of q 2 . 

On the contrary, in the annihilation NN — > 7* — > 
e + e~ the mass of virtual photon is equal to the total cm. 
NN energy. Its four-momentum squared is always greater 
than 4?w 2 . This allows to measure the nucleon form fac- 
tors in the time-like domain of q 2 , above the NN thresh- 
old. In this reaction, in order to study the form factor be- 



havior in a narrow domain near threshold, where non- 
trivial structures are predicted [9], one should have a beam 
of almost stopped antiprotons. This non-easy technical 
problem was solved at LEAR (10J. However, the under- 
threshold domain < q 2 < 4m 2 remains kinematically 
unreachable in this type of experiments. 
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Figure 2. The process pA — > (A — 1)7* (followed by 7* — > 

e+e-). 

One can penetrate in this domain of q 2 in the p anni- 
hilation on nuclei 

pA^ (A- 1) e + e~~ , 

see fig. 12 The symbol (A — 1) means not necessarily a 
nucleus but any system with the baryon number A — 1. 
Since extra energy of the antiproton can be absorbed by 
the (A — 1) system, the e + e~ pair may be emitted with 
very small invariant mass. Therefore the two-body reac- 
tion pA — ► (A — 1)7* is kinematically allowed for a very 
wide domain of invariant mass of the 7*, which starts 
with two times the electron mass, namely: 

4m 2 < q 2 < (y/sjX- M A _i) 2 . 

One can achieve near-threshold, under-threshold and even 
deep-under-threshold values of q 2 even for fast antipro- 
tons. This however does not mean that this reaction pro- 
vides us direct information about the nucleon form fac- 
tors. For the latter, we should be sure that the observed 
e + e _ pair (and nothing more) was created in the annihi- 
lation pp — > e + e~ on the proton in the nucleus, i.e., that 
the reaction mechanism is given by the diagram of fig. |3] 
or by a similar diagram where the p can rescatter before 
annihilation. 

At the same time, since the nucleons in the nucleus are 
off-mass-shell, the form factors entering the amplitude 
of fig. |3) are not precisely the same as found in the free 
pp annihilation. In general, the three-leg vertex, shown 
in fig. |U depends not only on the photon virtuality q 2 , 
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Figure 3. Amplitude of the reaction pA — > ( A — 1)7* in impulse 
approximation. 

but also on the nucleon ones p?, p\. F = F(p 2 , p 2 , q 2 ). In 
the case considered, the incident antiproton is on-ener- 
gy-shell: pp = m 2 , however the form factors depend on 
the proton virtuality p 2 . How the form factor F{pp = 
m 2 ,p 2 7^ m 2 ,q 2 ) v.s. q 2 differs from the free one F(p? = 
m 2 , p 2 = m 2 ,q 2 ) - this depends on the dynamics deter- 
mining its behavior v.s. the nucleon leg virtuality The 
nucleon form factors with off-shell nucleons were stud- 
ied in the papers 116111.71 . When only one nucleon is off- 
shell, then there are six form factors instead of two in the 
on-shell case. Generally, we can expect that the form fac- 
tor dependence v.s. p 2 , is much smoother than the q 2 de- 
pendence. The p 2 „ dependence can be determined by the 
nucleon self-energy corrections (i.e., by the structure of 
the nucleon), whereas the q 2 dependence in the time-like 
domain is governed by the pp interaction. The nucleon 
dynamics has a much larger energy scale than the nu- 
clear one. The typical off-shell variation found in the pa- 
pers II16II17| was from a few to 10 percent. We do not pre- 
tend to such an accuracy here. Therefore we neglect this 
effect in our calculation. We will come to this question 
later. In any case, both domains: q 2 < Am 2 , p 2 = m 2 and 
q 2 < Am 2 , p 2 < m 2 are totally unexplored experimen- 
tally and are interesting and intriguing. 



Q 2 



Figure 4. Nucleon electromagnetic vertex. 



We emphasize that though the form factor dependence 
on p^ can be weak, the nucleon off-mass-shell effect is 
very important for the kinematical possibility to reach 
the near- and under-threshold domain of q 2 with fast an- 
tiprotons. To produce the near- threshold e + e _ pairs in 
annihilation of a fast p on an on-mass-shell proton, the 
antiproton should meet in the nucleus a fast proton with 
parallel momentum. The probability of that, which we es- 
timated in the Glauber approach, is negligibly small [18] 
relative to the results presented below. However, if the 
effective mass p 2 , of the virtual proton is smaller than 
m 2 (that is just the case in a nucleus), then the near- and 
under-threshold e + e~ pairs can be produced in collisions 
with not so fast intra-nucleus nucleons. This effect con- 
siderably increases the cross section. To have an idea of 
the order of magnitude which one can expect for this 
cross section, we will calculate it in the impulse approx- 
imation. Numerical applications will be done for the low- 
est antiproton beam momentum foreseen in future projects. 
Namely, at the High Energy Storage Ring at FAIR-GSI 
this value is 1.5 GeV/ c. 

2.1 Cross section calculation 

At first, we consider the case of the deuteron target. If 
we know the amplitude of the reaction pd — > e + e~n: 
Mp^ e + e - n (to be calculated below), then the correspond- 
ing cross section is given by: 

^ cr pd->e+e~n = ' JT \^pd^e+e-n\ 2 
xS^(p p + p d - p e + - p e - - p n ) 

d 3 p e + d 3 p e - d 3 p n 
(2n) 3 2e e+ {2nf2e e - (2n) 5 2e„ { ' 

where I results from the flux factors: 
I = Fp£d = {{Pp ■ Pd) 2 ~ m 2 Mj} 1/2 = M d p pM « 2mp pM , 
PpXab is the incident p momentum in the lab. system. Here 
and below we imply the sum over the final spin projec- 
tions and average over the initial ones. 

We are interested in the distribution in the invariant 
mass M. of the final e + e~ system. To find it, for fixed 
value of Ai, we can integrate, in some limits, over the 
angles of the recoil neutron (determining the neutron re- 
coil momentum) and over the angles of the emitted e + e~ 
in their center of mass. This can be done using standard 
techniques of the phase volume transformations. Namely, 
we use the identity: 

/ S {i) {V e+e - - Ve + - p e -)d 4 V e+e -S(V 2 +e - - M 2 )dM 2 = 1, 

where V e + e - is the four-momentum of the e + e~ pair and 
.M 2 is its invariant mass squared. We insert this formula 
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in the cross section (0 and represent the three-body phase 
volume in 10 as: 



dc, 



n = /...^ (4) (Pp + 



pd— >e+e 



d 3 p e + d 3 p e - d 3 p n 



Ve+ ~ Ve 



Pn) 



2e e + 2e e - 2e n 
= J ...S^( Pp + p d -V e+e --p n ) 
xS i4 \V e+e - - p e+ - Pe-)d 4 V e+e -S(r^ +e - - M 2 )dM 2 
...dVi dV 2 dM 2 



d 3 p n d 3 p e+ d 3 p e - _ 



2e n 2e e + 2e e - 

where e n is the energy corresponding to the four-mom- 
entum p n and similarly for other energies, and we give 
explicitly only the factors resulting from the phase vol- 
ume. We denoted the phase volume of the e + e~ pair as: 



,-) 



d 3 p e - d 3 p e 



2e e - 2e e + 



AM 



dClp ~ — dn P 



(3) 



Here p* is the final electron (or positron) momentum in 
the cm. frame of the e + e~ pair and f\ is its solid angle 
in this frame. Neglecting the electron mass, we replaced: 

p* e « \M. 

The two-body phase volume dVi of the final nj* state 
reads: 



dV, 



J S (4 Hpp + p d -V e+e --p n ) 



d 3 V e+e - d 3 p n 
2e M 2e n 



7*n 



4,/s 



da 



pd 



(4) 



Here pl* n and n** n are in the cm. frame of the reaction 
(i.e., the cm. frame of pd) and Sp d is the corresponding 
cm. energy squared. 
In this way we obtain 



da, 



pd—tC 1 



2 8 7T 5 



mppjab 



\M 



pd^e+e 



- n \ 2 dV l dV 1 dM 2 . 

(5) 

This form of the cross section will provide us the invari- 



ant mass distribution dcrp d ^, e + e - n / dM. 

The formula ((5} is general and does not assume any 
particular mechanism of the reaction. In the impulse ap- 
proximation, when the mechanism is given by the dia- 
gram of fig.|3j the total amplitude squared \Mp^ e + e - n \ 2 
is proportional to the annihilation amplitude squared 
| Mpp^ e + e - 1 2 and to the square of the deuteron wave func- 
tion \ip\ 2 . The direct calculation of the amplitude corre- 
sponding to the diagram of fig. Ogives the proportional- 
ity coefficient: 



|M 



pd—>e + e n I 



4m|Mp p _ >e+e - 



M 2 , 



(6) 



and \ ip\ 2 is normalized as: 



m 



d 3 k 



1. 



(2n) 3 e k 

Substituting ^ into ©, we find the cross section: 



(7) 



da, 



Vp,\ab 



\-^-pp—*e + e 



\ 2 ^ dn e 



x|4ffpzrfn; tn dM 2 . 



We integrate over dCl e (in the e + e cm. system) and fi- 
nally obtain: 



da. 



pd- 



dM 



' pp^e+e 



(M)t](M), 



(8) 



where t](M) is the distribution (given by eq. (flTb below) 
of the e + e~ invariant mass M due to the fact that the re- 
action occurs on the proton bound in the nucleus, and 
Cp p ^e+ e -(M) is the cross section of the pp — > e + e~ an- 
nihilation at the total energy M . The latter is determined 
by the amplitude Mp^ e + e - , corresponding to annihila- 
tion via the s-channel photon 7* . The calculation of this 
amplitude is standard. Though the target proton is off- 
mass-shell, we do not consider the off-mass-shell effects 
in this amplitude, i.e. we assume a free proton. Then the 
amplitude squared reads: 



|M 



pp- 



16a 2 n 2 
M 2 

4l\r. J2 



[4m 2 |G £ | 2 sin 2 t , + 



M l \G M \ 2 cos 2 0, 



e+\> 



(9) 



where 6 e + is the emission angle of the e + or the e~ in 
the 7* cm. frame and G £ , Gm are the proton electric and 
magnetic timelike form factors. To estimate the nuclear 
effect, from now on we omit the form factors, i.e. we put 
I Ge I = \Gm\ — 1/ or = 1/ IF2I = as for a pointlike 
nucleon. Then the pp — > e + e~ cross section obtains the 
form: 



PP- 



647r 2 mpp i ab M 

2a 2 n(2m 2 + M 2 ) 
3M 2 mpp Mb 



\M 



PP- 



L dC\ e 



(10) 



Though crpp^ e + e - depends on M, the main (nuclear) 
effect is determined by the factor rj (M ) : 



V (M) 



mi 



(2ny 



spd 



\ip(k)\ 2 dz- 



(11) 



Here z = cos 8, where 8 is the angle, in the cm. frame 
of the reaction, between the initial deuteron momentum 
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p* d and the final neutron momentum p*. The integration 
over dz in f](M.) results from dO* tfl in the phase volume 
dV2, eq. I0J. The argument of the wave function k depends 
on z. This explicit dependence is given in the next section. 
The distribution rj(AA) is normalized to 1: 



i]{M)dM = 1. 



(12) 



dilation of a distribution does not change its normaliza- 
tion: t](A4) remains normalized to 1. 

For a test of calculation we can find from 10 (or ex- 
tract from (flOt ) the amplitude squared \Mp f) _ >e + e - | 2 , av- 
eraged over the angle 6 e + : 



64n 2 ci 2 (M 2 + 2m 2 ) 
3M 2 



To prove eq. CG), we calculate this integral explicitly: Due to time invariance, it is equal to | M e+e -^ pp \ 2 . Then 



tj(M)dM 



(2n) 



*pd 



Since 



(s pd -M 2 + m 2 )/{2, 



we have: M dM 
this way we find: 



Sprfde*. Then: de*, 



In 



(13) 

Pn d Pn J n 



t]{M)dM = / \ip(k)\ 



= J \m\- 



> mp* n 2 dp* n 
(2tt) 3 4 
, md 3 k 



dCl* 



The latter equality is just the normalization condition (0 
of the wave function. Here the neutron momentum p* is 
defined in the cm. frame of the total reaction, whereas 
the momentum k is defined in the cm. frame of the np 
system. So, p* n and k are the momentum of the same parti- 
cle, but in different frames. We replaced mp* 2 dp* n dd* / e p * 

by ^^M, since this integration volume is a relativistic in- 
variant. 

The total cross section is obtained by integrating I0 
in the finite limits J^A-rnvn 

< M<M max* where M.fnin — 
2m e , Mmax = ^/spd ~ m - Neglecting the electron mass, 
we can put M. m {„ = 0. For p*/^ = 1.5 GeV/c the value 
M max is high enough and provides the normalization 
condition (121 1 with very high accuracy. 

To emphasize more distinctly the effect of the nuclear 
target, we can represent the cross section (10)| of the an- 



nihilation pp 
©: 



e + e on a free proton similarly to eq. 



f'P- 



dM 



a, 



pp- 



e-HM 



(14) 



where cr. 



PP- 



is defined in i fTOl l, s 



PP 



(Pp + Pp) 2 



2m 2 + 2m s Jp 2 M - 

The fact that in the annihilation on a free proton the 
mass of the final e + e~ pair is fixed is reflected in iQj} in 
the presence of the delta-function. Comparing this for- 
mula with lO, we see that the effect of the nuclear tar- 
get results in a dilation of the infinitely sharp distribution 
S(Ai — y/Spp) in a distribution of finite width rj(A4). The 



we can identify m with the muon mass and find the total 



cross section of the reaction e + e 



c + 



*}^P" 



6AnE 



*3 



(p* and E* are the final muon momentum and energy). It 
coincides with a well-known result and, at high e + e~ en- 



ergy, turns into a e + e 



Ana 1 



given in textbooks. 



This confirms the correctness of numerous factors appear- 
ing in our formulas. 



2.2 Analysis and numerical calculations 

So far we have not defined the argument k in the wave 
function ip(k) determining the distribution i jTTV The inte- 
gral over dz incorporates an interval of k, which, as we 
will see, for pp = 1500 MeV/ c and M. = 2m starts with 
the minimal value fcmin « 360 MeV/ c. For these values 
the relativistic effects become non-negligible, therefore dif- 
ferent relativistic approaches result, in principle, in dif- 
ferent expressions for k. If we take as the argument the 
spectator momentum (the momentum of the neutron) in 
the rest frame of the deuteron, then it is expressed as: 
k 2 = e 2 — m 2 , where the neutron lab. energy can be rep- 
resented as e„ =(pd' Pn)/Md, so that when p d = (M^, 0) 
we obtain identity. Introducing the invariant t = (pd — 



p„) 2 , and since p d ■ p n = (M; 



t)/2, we get: 



(M 2 + m 2 - t) 2 2 

k = — = m . 

AM 2 



In the cm. frame of the total reaction pd 
variable t is expressed as: 



(15) 



e+e n, the 



t = (Pa 



Mi 



2(e* d e* n 



zp*dP*n)+m 2 (16) 



where e* d , e* and p* d , p* n are the cm. energies and momenta 
of the deuteron and the neutron, z = cos 8 and 6, as al- 
ready mentioned, is the angle between p *, and p* n . The en- 
ergy e* n is given by eq. l fl3)l , whereas e* d is obtained from 
| [T3l by the replacement of Ai 2 by M 2 . After these preci- 
sions, the formula l fl5|l together with | [T6)| completely de- 
termines the argument of the wave function. 
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It should also be emphasized that the proton is off- 
mass-shell. Its off-shell mass squared has the value (m* ) 2 = 
p 2 = f, which depends on kinematics. 

The minimal value /c m in is achieved at z = 1, i.e. when 
the neutron recoils antiparallel to the incident antiproton. 
Using eqs.dl5t and fT6t , one finds a value that is remark- 
ably small: pa 360 MeV/ c, corresponding to an off- 
shell mass m* pa 800 MeV (pa 0.85m) for the initial proton. 

We emphasize the meaning of the integration over z 
in eq. {TT)l . The incident antiproton with momentum pp in 
the pd cm. frame meets in the deuteron the proton with 
momentum p*. In the impulse approximation (diagram 
of fig. |3}, the latter is determined by the difference of the 
deuteron and the neutron momenta: p* — p\ — p* v Since 
eq. CO} gives the value of rj (Ai) at fixed Ai, for all these 
momenta and the relative angles between them and for 
the virtual proton mass squared (m*) 2 = t, the invariant 

2 is just equal to the given in- 



of eq.lfTBll, reducing the minimal value fc n 
only. 



by 10 MeV/c 



energy ^fs^ = ^(p p + p p ) 

variant mass Ai of the virtual 7* and of the final e + e~ 
pair. In this way, we obtain the contribution of the events 
with given time-like value of q 2 = Ai 2 in the cross sec- 
tion. When the angle 9 between p* and p* d (and z = cos 9) 
varies, the modulus of the proton momentum, its direc- 
tion and the off-shell mass m* vary correspondingly to 
provide the same fixed value of Spp = Ai 2 = q 2 . The 
integral over z in fTTT l incorporates all these events, with 
different p*, m* and with the same q 2 . 

At first glance, the small near-threshold pp cm. en- 
ergy Ai pa 2m in the collision of a fast p (pp = 1500 
MeV/c) is achieved, when the antiproton meets in the 
deuteron a fast proton having the same momentum as 
the p, in modulus and direction. The protons with such 
a high momentum are very seldom in deuteron. For this 
mechanism, the cross section would be very small. How- 
ever, the near-threshold value of Ai is obtained in other 
kinematics. As we mentioned, the proton momenta k in 
the deuteron wave function ip(k) in eq. fTl) start with 
~ Knin ~ 360 MeV/c only (that corresponds to z pa 1). 
The main reason which allows to obtain in this collision 
the value Ai ~ 2m is the off-shellness of the proton: 
m* < 0.85/7? instead of m* = m. This 15% decrease rel- 
ative to the free proton mass is enough to obtain the in- 
variant pp mass Ai pa 2m, when one has the two parallel 
momenta: 1500 MeV/c for p and 360 MeV/ c for p. 

We also notice that for other values of z, i.e. non-colinear 
p and p momenta, the proton momentum k needed to 
produce Ai pa 2m is larger than k m [ n . In a correlated way, 
this proton is also further away from the mass-shell. 

The calculation in the framework of another relativis- 
tic approach - light-front dynamics fl9ll - gives another 
formula for k which numerically is very close to the one 




1800 1900 2000 2100 2200 2300 

M (MeV) 



Figure 5. The cross section 



Ac, 



1W 



- of the reaction pd 



e+e~n v.s. M, in the interval: 1750 MeV < M < 2350 MeV, 
calculated for a pointlike proton. The arrow indicates the pp 
threshold. 

The cross section d<7p d ^ e+e - n / dM, eq. ©, has been 
calculated for an antiproton of momentum pp = 1500 
MeV/c on a deuteron nucleus at rest, with the deuteron 
wave function Il20l , incorporating two components cor- 
responding to S- and D-waves. The result is shown in 
fig. [5] The maximum of the cross section is at Ai = 2257 
MeV, that corresponds to the p interacting with a proton 
at rest (and on-shell). The cross section integrated over 
Ai is equal to 43 nb. We remind that these calculations 
do not take into account the proton form factor. Its influ- 
ence will be estimated below. The numerical integral over 
Ai of the function rj(Ai), eq. ((Til l, is ~ 1, in accordance 
with the normalization condition (Tl2l l. 

The pp threshold value Ai = 1880 MeV is on the tail 
of the distribution, far from the maximum. Relative to the 
maximum, the cross section at threshold decreases ap- 
proximately by a factor 600. The numerical value at the 
threshold is: 



da 



dAi 



M=2m 



pb 
MeV 



(17) 



In fig.[6]this cross section is shown in the near-threshold 
interval 1830 MeV < Ai < 1930 MeV. The integral over 
Ai in a bin of width 100 MeV centered on the threshold 



is: 



1930 MeV da 



pd- 



dAi pa 100 pb . 



1830 MeV dAi 

These estimations take into account the suppression 
resulting from the momentum distribution in deuteron. 
However, they do not incorporate the form factors of the 
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Figure 6. The same as in fig.|5] but near the pp threshold, in the 
interval 1830 MeV <M< 1930 MeV. 



nucleon. To incorporate them in a simplified way, one can 
consider an effective form factor |F| which depends on 
Ai, and include it in the integral: 




M (MeV) 



Figure 7. The cross section — r> 'dM " mcru ding an effective 
form factor (see text) is drawn as a dashed curve above the pp 
threshold. The solid curve is the same as in fig. [5] (but now in 
logarithmic scale). 

2.3 Annihilation on heavier nuclei 



>pd- 



j a vv- 



(M) t](M) \F(M)\ 2 dM (18) 



where <Tpp^ e + e - (Ai ) is the cross section for pointlike nu- 
cleons given in eq. (p~0b . To have an estimate of this inte- 
gral, we have taken the effective proton form factor mea- 
sured in ref. IfTTI . By doing this, we neglect all off-shell 
effects. We interpolate \F(Ai) | linearly between the mea- 
sured values, and we limit the integral to the region Ai > 
2m. In this way we obtain Cp^e+e-H — 1 n b, which is 
comparable to the total cross section ap„^ e + e - on a free 
proton at M = 2257 MeV.0 

Figure [7] shows the integrand of eq.([18l as a function 
of Ai, using the above choice of form factor. One notices 
the rapid rise of the cross section as Ai approaches the pp 
threshold. Lastly, we point out that at threshold, our dif- 
ferential cross section dap^ne+e- of 1 pb/MeV (eq.llT7Tl) 
is not suppressed by any factor, since there the form fac- 
tor |F| seems to be close to 1 experimentally I21II22I . Be- 
low this threshold one may expect a form factor effect 
larger than one, but this is not known and this is the goal 
of the proposed study. 



This cross section u. 



pp—*e + e~ 



1 nb at M 



obtained from the measured cross section cr e+e -^pp 



2257 MeV is 
at the same 

cm. energy, which we estimate to be ~ 318 pb, from interpo- 
lation in the data of ref. I'll |. This value is then multiplied by 
the factor M 2 / (M 2 — Am 2 ) w 3.27 to obtain the cross section 
in the inverse channel pp — > e + e~ at the same cm. energy. 



For A > 2, we should take into account the possibility of 
excitation and breakup of the final nucleus A — 1 in the 
process pA — > (A — 1)7*. The result contains the sum 
over the final energies of the residual nucleus and the in- 
tegral over a continuous spectrum. That is, the function 
I if (it) | 2 ineq. flU is replaced by the integral fj^™ S(E, k)dE, 
where S(E,k) is the nucleus spectral function giving the 
probability to find in the final state the nucleon with the 
relative momentum k and the residual nucleus with en- 
ergy E. For high incident energy we can replace the upper 
limit by infinity. Then we obtain: 



S(E,k)dE = n(k), 



where n(k) is the momentum distribution in the nucleus. 

To estimate the cross section on heavy nuclei, we will 
still use eqs. (|8), i fTTll but with the two following changes. 

(i) We replace the deuteron momentum distribution by 
the nuclear one. Since the deuteron wave function is nor- 
malized by (0, whereas n(k) is usually normalized as 
J n(k)d 5 k = 1, we replace ip 2 (k) in ((TlT lby (2n) 5 ei c n(k)/m. 

(ii) We multiply iQj} by the number of protons Z. 

The numerical calculations were carried out for the 
12 C, 56 Fe and 197 Au nuclei with the nuclear momentum 
distributions found in the papers 1 23 , 24 ] . Near threshold, 
i.e. at M. = 1880 MeV, for all three nuclei we obtain very 
close results given by: 



d~M 



6.5 Z 



pb 
MeV 



(19) 
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Multiplying by the charge Z (Z( 12 C) = 6, Z( 56 Fe) = 26, 
Z( 197 Au) = 79) and integrating (19) over a 1 MeV inter- 
val near M = 1880 MeV, we get: 

cr( 12 C) = 39 pb, cr( 56 Fe) = 0.17 nb, 0-( 197 Au) = 0.5 nb. 

These results were found without taking into account 
the absorption of p in nucleus before electromagnetic an- 
nihilation. This absorption was estimated in Glauber ap- 
proach (it is applicable since now we do not need high 
nucleon momenta). It reduces these cross sections by only 
a factor 2. We recall again that these results are obtained 
for structureless nucleons. 



and the neutron ones. One can expect that the direct anni- 
hilation pp — > e + e~ dominates over the mechanism with 
preliminary transition pp — > fin. In this case, the contri- 
bution of the proton form factor dominates. 

We emphasize that in any case, whatever the inter- 
mediate steps are in process {I}, the e + e~ pair of the final 
state must come necessarily from the baryon-antibaryon 
electromagnetic annihilation, pp or nn, because there is 
only one neutron left at the end. It cannot come from an- 
other process, even if there are complicated intermedi- 
ate steps, like rescattering, etc. Therefore this e + e~ pair 
is a direct and very little distorted probe of the baryon- 
antibaryon electromagnetic annihilation vertex. 



2.4 Beyond the impulse approximation 

Without carrying out any calculation, we discuss in this 
section other possible mechanisms for the process pd — ► 
e + e~n. One of them is the initial state interaction, which 
includes rescattering (not only elastic) of the initial p in 
the target nucleus. In the rescattering, the incident p looses 
energy and therefore the proton momentum needed to 
form the invariant mass A4 ~ 2m becomes smaller. The 
probability to find such a proton in deuteron is higher. 
Therefore initial state interaction increases the cross sec- 
tion. 

Though, if the rescattering is inelastic (with pion cre- 
ation), it results in the reaction pd — > e + e~Nn. In this 
reaction, the e + e~ pair is still produced in the pp electro- 
magnetic annihilation, giving an information about the 
proton timelike form factors. Therefore this process is also 
interesting in itself, although it is not in the scope of the 
paper. 

Another mechanism of the reaction pd — ► e + e~n may 
include also the transition pp — ► fin with the subsequent 
annihilation fin — > e + e~ . This transition can take place, 
for example, via meson exchange: 

pp — > {nn~)p — > n(n~p) — > fin, 
pp — ► p(n + n) — ► (pn + )n — > nn. 

If the annihilation fin — > e + e~ occurs on the same neu- 
tron which was created in the reaction pp — > nn, this does 
not give anything new, since this is simply a particular 
contribution pp —> fin — ► e + e~ in the initial state interac- 
tion incorporated in the full amplitude pp — ► e + e~ , i.e., 
in the timelike proton form factor. 

On the contrary, when the annihilation nn — > 7* — > 
e + e~ occurs on another neutron (the neutron from deu- 
teron), then the reaction pd — > e + e~n gives information 
about the timelike neutron form factor. So, what is mea- 
sured is the sum of two timelike form factors: the proton 



3 Experimental aspects 



The aim of this section is to investigate the feasibility of 
process (TJ in future experiments using antiproton beams. 
The study is made in the case of the PANDA experiment, 
which will use a state-of-the-art internal target detector at 
the High Energy Storage Ring (HESR) at FAIR-GSI. This 
study is only a first step, exploratory and rather qualita- 
tive. It remains to be pursued in more details in the fu- 
ture. 

Regarding the target nucleus, we restrict ourselves to 
the deuteron. It is the only nucleus yielding a simple, 
three-body final state (e + e~n), for which an experimen- 
tal strategy can be defined. For heavier target nuclei A, 
the break-up channels will dominate over the three-body 
final state e + e~ (A-l), and no simple identification strat- 
egy can be defined. Indeed, to be able to use kinemati- 
cal constraints in order to select the desired reaction, one 
needs to detect (at least) all final particles but one, a task 
which becomes more and more difficult as the number 
of nuclear fragments increases. As a side remark, we also 
note that the HESR luminosity in (pA) decreases with the 
atomic charge Z of the target nucleus [25J, in a way that 
roughly compensates the increase of cross section with Z 
reported in sect. 12.31 

New technologies for antiproton accelerators, targets 
and detectors provide a gain in luminosity and allow to 
access processes with very low cross sections, like the one 
presented here in the near-threshold region. In this sec- 
tion we first examine reaction (the signal) and then the 
competing background processes. The presented mate- 
rial largely relies on the design performances of the PANDA 
detector [25 .26] and the study of the electromagnetic pro- 
cess pp^e + e"ona free proton at rest [25.26.27j. 
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3.1 Beam and target operating conditions 

Considering the process p + Pbound ~~ > 1* on an off-shell 
proton bound in the nucleus, the lower the antiproton 
beam momentum, the easier it is to reach a photon vir- 
tuality Q 2 < Am 2 (Q 2 = q 2 ). Ideal conditions would be 
antiprotons of very low momentum. Here we consider 
the nominal lowest beam momentum that is foreseen in 
the HESR, i.e. Ppj a b = 1.5 GeV/c. As seen in sect. 12.21 in 
these conditions only a very small fraction of the target 
protons in the deuteron will have enough off-shellness in 
order to produce a pp system of invariant mass squared 
< 4m 2 . In the impulse approximation, i.e. assuming that 
the neutron is a spectator in the pd — > e + e~n process, 
these off-shell protons are found in the upper tail of the 
momentum distribution n(k) of the deuteron, at k = 360 
MeV/c and above. 

In the PANDA experiment, nuclear targets will be used 
to study hadrons properties in the nuclear medium. For 
deuterium, an internal target of cluster-jet type or pel- 
let type will be placed in the beam, and an effective tar- 
get thickness of 3.6 ■ 10 15 atoms/cm 2 is considered. At 
Pp,lab = 1-5 GeV/c, beam losses limit the number of an- 
tiprotons to 10 11 per cycle, yielding a maximal value of 
5 • 10 31 cm .s for the cycle average luminosity C in 
(pd) collisions 11251 . 

3.2 Count rate estimate for the process pd — ► e + e~n 

With the above operating conditions (ppj a b> Q an d me 
cross section a calculated in sect. 12.21 one can compute 
the number of events from process CU in an ideal detec- 
tor, covering the Arc solid angle with 100 % efficiency; it 
is N = C ■ <j . For one month effective beamtime (i.e. cor- 
responding to an integrated luminosity equal to £ x 2.6 • 
10 6 s), this yields a large number of events integrated over 
the full range of dilepton invariant mass: N ~ 1.3 x 10 5 
(this number is deduced from the cross section CTpp^ e + e - ~ 
1 nb on a free proton at rest). However, in the special 
region of interest near the NN threshold, the cross sec- 
tion da/dAi is much reduced, of the order of lpb/MeV 
(see sect. 12.2b , yielding a rate of 130 events per month in 
a 1 MeV bin of dilepton mass. This number of events is 
small, but other rare processes considered in the PANDA 
physics program have comparable expected rates. 

On the one hand, the global detector efficiency and 
the necessary experimental cuts will reduce this amount 
of good events, with a reduction factor that could reach 
50%. On the other hand, the cross section calculated in 
sect. 12.21 for process {1} is probably a lower limit in the 
subthreshold region, since there the proton form factor 



can be larger than one and enhance the cross section. There- 
fore, all in all, within one month effective beamtime one 
could obtain an experimental spectrum of dilepton mass 
Ai containing ~ a hundred signal events per MeV bin 
near the threshold (Ai ~ 2m). This would allow a first 
insight into possible structures in this region, including 
the totally unexplored subthreshold region. 

3.3 Experimental signature of the exclusive channel 
pd — > e + e~n and related background processes 

The selection of the exclusive reaction (Q} is based on the 
detection of the lepton pair. The experimental strategy re- 
lies on three main characteristics of the detector: 1) its res- 
olution, 2) its hermeticity and 3) its particle identification 
capability. 

Regarding the first point, the main variable is the miss- 
ing mass Mx, i.e. the invariant mass of the missing sys- 
tem X in the reaction pd — > e + e~ X. It is defined using the 
four-momentum vectors of the initial and detected parti- 
cles: 

M X = (Pp + Pd- Pe+ ~ Pe- f ■ 

Due to the many processes producing a lepton pair inclu- 
sively, this missing mass will have a wide distribution; 
for the events of reaction IQ} a peak must be searched at 
the neutron mass, which represents the physical lower 
bound of the Mx spectrum. Therefore the experimental 
resolution in this missing mass is a crucial parameter (see 
sect. 111). 

The background due to inclusive lepton pair produc- 
tion is quite large at the raw level. True e + e~ pairs orig- 
inate mostly from Dalitz decay (or direct l + l~ decay) of 
mesons, and also from real photon conversion in the tar- 
get. These types of processes will give dileptons mostly at 
low invariant mass (Ai < 1 GeV), while we are looking 
for high-mass ones (Ai > 1.8 GeV). Uncorrelated pairs 
coming from combinatorial background may form a high 
invariant mass; however they can be subtracted using the 
like-sign pairs (also combinatorial) as in heavy-ion dilep- 
ton experiments. All these processes of inclusive lepton 
pair production correspond to missing masses Mx larger 
than one neutron mass, and a large fraction of them can 
be eliminated by a proper cut in Mx- Also, these pro- 
cesses create more than three particles in the final state, 
and a condition on the observed particle multiplicity should 
be efficient to reject a large fraction of this background. 

This brings us to the second important aspect, of de- 
tector hermeticity. The PANDA detector will cover al- 
most the An solid angle, by combining a target spectrom- 
eter for large polar angles and a forward spectrometer 
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for small angles. The good hermeticity is essential in the 
physics program in order to detect a complete spectrum 
of final states. In reaction (I} the neutron will go unde- 
tected most of the times; actually in our kinematics of 
interest near threshold, this neutron is emitted at very 
backward angle w.r.t. the beam and will not hit any detec- 
tor. With respect to hermeticity potential background to 
reaction ([]} consists in reactions in which, in addition to 
the e + e~ pair, one or more light particles (pions, photons) 
are emitted and escape detection. They will contribute to 
form a missing system X (as defined above) of low bary- 
onic mass, typically Mx — m + m n , polluting the region 
of interest. The simplest cases of such background reac- 
tions are: 

pd e + e~n n° (from pp bound -> 7* n°) 
pd^ e+e~pn~ (from pn bound -> 7* n~) 
on the bound proton or the bound neutron in the deuteron. 
The final tzN state can be non-resonant or resonant (e.g. 
the A(1232)). 

A third important aspect is the particle identification 
capability of the detector (PID), for charged and neutral 
particles. In PANDA, a global estimator ||25l combines 
the information of the various subdetectors and will pro- 
vide charged particle identification in the momentum rangi 
from 200 MeV/c to 10 GeV/c. In our case the most im- 
portant task is to discriminate between e and n . In- 
deed, the main competing background to reaction |TT) is 
of the type pd — > n + n~n or pd — ► tz + tt~X where the 
X system has a low baryonic mass, and the n + n~ pair is 
misidentified as an e + e~ pair. 

The case of p p — > e + e~ on a free proton at rest has 
been studied in detail E7I in view of future measure- 
ments of the timelike proton form factors in PANDA. 
The exclusive process pp — > n + 7Z~ is the dominant back- 
ground to this measurement, with a cross section ~ 10 6 
larger than the pp — ► e + e~ cross section in the HESR en- 
ergy range. Simulation studies have been performed in 
order to estimate precisely the fraction of pp — > n + 7T~ 
events that can be eliminated by PID cuts. The most effi- 
cient elements for tz^ /e discrimination are the electro- 
magnetic calorimeters and the dE/dx measurements in 
the tracking systems. Applying the global PID cut at the 
"very tight" levelH yields an efficiency to single electrons 
of ~ 90% with a contamination rate of pions smaller than 
10 -3 , for particle momenta above 1 GeV/c 1251 . As a re- 
sult, a rejection factor of 10~ 7 is reached on the pion pair 
of the process pp — > n + n~ . This background is further 
reduced by a factor 50-100 when performing a two-body 
kinematical fit. In this way one obtains a more than 99% 

2 corresponding to a global identification probability greater 
than 99.8% for an e + or an e~ . 



clean sample of p p — > e + e~ events, with a 20-30% global 
efficiency [27\. 

These results cannot be extrapolated directly to the 
case of reaction (T), for which a detailed simulation would 
be required. Nonetheless, the strategy is based on the same 
elements. In a simplified scheme, the use of the "very 
tight" PID cut allows to reduce the p d — ► n + n~ n hadro- 
nic background by a factor ~ 10 6 (10 3 for each charged 
pion). After that, one is left with a background- to-signal 
ratio of about 1:1. Kinematical constraints are less effi- 
cient than in the case of pp — > e + e~ since we have a 
three-body final state. Here only one constraint can be 
imposed, on the missing mass. Therefore the remaining 
background from hadronic processes cannot be eliminated 
by a cut, it has to be calculated and subtracted. It is in 
principle possible, provided that the cross section for these 
processes is measured (in PANDA) and the PID performan- 
ces of the detector are known with accuracy. Other re- 
actions due to strong pn annihilation are also potential 
sources of background. For example the process pn — > 
7r~7r° followed by rfi Dalitz decay, creates a n~e + pair 
that can be misidentified as an e~e + pair; etc. 

In summary, the simulations already performed for 
the channel pp — > e + e~ serve as a benchmark study to 
explore the possibility to measure the process p d — > e + e~n 
in PANDA. The competing background to this reaction is 
certainly more important than in the case of p p — ► e + e~ 
on a free proton at rest, and more difficult to reject. We 
have given elements for a strategy, but the full experi- 
mental feasibility remains to be proven. A further study 
of this point will require complex simulations which are 
beyond the scope of this paper. 

3.4 First-order simulation of pd — ► e + e ~n 

A simple Monte-Carlo simulation of reaction (Q}, of pure 
phase-space type, was performed. The detector resolu- 
tion is implemented in order to estimate some important 
parameters, such as the achievable resolution in missing 
mass (c"m x ) an d in dilepton mass (<7>f). 

Kinematics are defined by an incoming antiproton of 

1.5 GeV/c momentum onto a deuteron nucleus at rest. 
The two-body reaction p d — > ny* is generated by sam- 
pling in cos8j* cm (the polar angle) and in the azimuthal 
angle (p-y* C m m the pd center-of-mass. The photon vir- 
tuality Q 2 = (Ai) 2 is sampled between Q 2 njn = 4m 2 . 
and Qmax = (^/Sp^ — ni) 2 . Then the decay 7* — ► e + e" is 
generated in the 7* center-of-mass by sampling in cos 6 e - 
and <p e - . Finally all the produced particles are transformed 
back to the laboratory frame. The five samplings are made 
in uniform distributions, therefore the event weight is not 
realistic for this first study. 
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Figure 8. Simulation of the reaction pd — > e + e~n at antiproton 
beam momentum Ppi a j, = 1.5 GeV/c: electron energy versus 
positron energy in the Lab frame. The solid line corresponds to 
our threshold kinematics of interest (see text). The dashed line 
corresponds to the reaction pp-te + e"ona free proton at rest, 
at the same value of Ppj a b- 

The detector resolution is implemented on the mo- 
mentum and angles of the electron and positron produced 
at the vertex, by sampling a measurement error in a Gaus- 
sian distribution, independently for (p, Lab of each 
particle. For the relative momentum resolution (a p / p) we 
take a uniform value of 1.5% in r.m.s., which should be 
conservative given the rather low momenta involved (p < 
2.2 GeV/c). The angular resolution is taken to be eg = 
dtp = 3mr in r.m.s. 

The resulting phase space of the lepton pair is shown 
in fig. [8] for the energies and in fig. [9] for the polar an- 
gles. When compared to the case of the reaction p p — > 
e + e~~ on a free proton at rest (dashed curves), now most 
of the correlations between the two leptons are lost and 
the phase space is much more open. The special thresh- 
old kinematics of interest in reaction (l} is defined by: i) 
a photon virtuality corresponding to the NN threshold, 
i.e. A4 = 2m , and ii) the virtual photon emitted at for- 
ward angle in the lab (cf. sect. 12.21 this corresponds to 
z = 1 or k = k m j n in the deuteron momentum distribu- 
tion n(k)). For this particular kinematics, represented by 
a solid curve on Figs. [8] and |9j lepton-to-lepton correla- 
tions are re-established, but they differ from the case of a 
free proton at rest, because the target proton is now off- 
shell in the deuteron. 




50 100 150 



e Lab (e+) (deg) 

Figure 9. Simulation of the reaction pd — > e + e~n at ppj a i, = 
1.5 GeV/c: polar angle of the electron versus polar angle of the 
positron in the Lab frame. The solid curve corresponds to our 
threshold kinematics of interest (see text). The dashed curve 
corresponds to the reaction p p — > e + e~ on a free proton at rest, 
at the same value of Ppjat- 



The off-shellness of the initial proton in process Q} 
has been defined in sect. 12.21 In the impulse approxima- 
tion, the off -shell mass squared of this proton is equal to 
t = (pj — p n ) 2 , where pj and p n are the four-momentum 
vectors of the initial deuteron and final neutron. In this 
approximation, the proton off-shellness is an experimen- 
tally accessible quantity since the neutron is "measured" 
as the missing particle. Figure [10] shows the t variable 
as a function of the photon virtuality, in the full phase 
space. In the threshold kinematics of interest, f is equal 
to 0.63 GeV 2 (to be compared with the free proton case, 
t = M 2 p = 0.88 GeV 2 ). 

Finally, fig.lTTlshows the experimental resolution that 
can be expected on the most important variables: the miss- 
ing mass Mx and the dilepton invariant mass Ai . These 
resolutions are quite uniform within the reaction phase 
space, apart from a slight Q 2 -dependence. The width of 
the distributions in fig. [TT] is dominated by the momen- 
tum resolution of the detected e + e~ pair. The obtained 
resolution in missing mass is rather good, with an r.m.s. 
of 46 MeV, or a full- width at half -maximum of 106 MeV 
This value is smaller than one pion mass, therefore by 
applying a cut around the neutron mass in the Mx spec- 
trum one should be able to separate - at least partly - the 
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Figure 10. Simulation of the reaction pd — > e + e~n at Ppj a b — 
1.5 GeV /c: the variable t (see text) versus the photon virtual- 
ly Q 2 — (M) 2 . The circled star corresponds to our threshold 
kinematics of interest, where t = 0.63 GeV 2 . 



Figure 11. Simulation of the reaction pd — > e + e~n at ppj a i, = 
1.5 GeV/c: resolution obtained for the missing mass Mx (upper 
plot) and for the dilepton invariant mass M. (= M e + e -, lower 
plot). 



missing systems X of the type (one nucleon + one pion), 
which are the closest in the spectrum. 

The resolution in dilepton invariant mass (r.m.s. = 17 
MeV) is also reasonably good, and well suited to evi- 
dence some structures due to possible baryonium bound 
states near the pp threshold. 

4 Other aspects 

In this paper we have concentrated mostly on the pos- 
sibility to use process Q} to access the pp threshold and 
under-threshold region. However, we would like to stress 
that this process is also interesting for other purposes. 

First, it provides a full range of e + e~ invariant mass 
M. with one single antiproton beam energy. In this sense, 
the role of the nucleus is a bit similar to the role of initial 
state radiation in the inverse reaction e + e~ — > pp7/sR.As 
a consequence, any clean experimental data on process 
Q} at any beam energy (e.g. HESR momenta higher than 
1.5 GeV/ c), have the potential to yield valuable informa- 
tion on the (half off-shell) nucleon timelike form factors, 
in a large Q 2 -range. Of course, in this perspective more 
theoretical work is needed to test the validity of the im- 
pulse approximation, and to better define the off-shell ef- 
fects in the NNj* vertex. 

Second, one should note that the antiproton momen- 
tum 1.5 GeV/ c just corresponds to the threshold value of 



creation of the AA pair on a free proton. Therefore the 
virtual creation of the AA pair in reaction (l} is not sup- 
pressed by the nucleon momentum distribution in deuteron 
and contributes just in the domain of the peak of fig. |5J 
that allows one to study the AA threshold region with 
good statistics. In the AA system, the quasi-nuclear states 
were predicted in [28 [ and, similarly to the NN quasi- 
nuclear states, they should manifest themselves as irreg- 
ularities in the cross section. The contribution of the chan- 
nel pp — > A A — > e + e~ in the total cross section pp — > 
e + e~ (the latter equals 1 nb, see sect. 12.21 above) was es- 
timated in l29l as 0.1 nb, i.e. 10% of the total cross sec- 
tion. Therefore we expect that the structures caused by 
the channel pp — > AA — > e + e~ can be observed in pro- 
cess (1) in the region of mass A4 near the AA threshold. 

5 Conclusion 

We have studied the reaction pA — > (A — (followed 
by 7* — > e + e~). This process gives access to the pp an- 
nihilation pp — > 7* at invariant masses v^pp which are 
below the physical threshold of 2m, due to the proton 
off-shellness in the nucleus. In this way a possibility ex- 
ists to access the proton timelike form factors in the near- 
threshold and the totally unexplored under-threshold re- 
gion, where NN bound states are predicted. 
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The differential cross section da/dAi has been cal- 
culated as a function of the dilepton invariant mass Ai, 
for an incident antiproton of 1.5 GeV/c momentum on 
a deuteron target (and heavier nuclei). The calculation is 
done in the impulse approximation; then the distribution 
of Ai is obtained from the deuteron wave function ip(k). 
We find that the pp threshold (Ai = 2m) is reached for a 
minimal proton momentum ^^=360 MeV/c in the nu- 
cleus, and at this point the cross section is about 1 pb/MeV. 
The calculation does not include the form factor effect, 
which should come as an extra factor. 

Experimental aspects have been investigated in the 
case of a deuteron target, i.e. for the three-body process 
pd — > e + e~n. We have taken the conditions of the PANDA 
project at FAIR-GSI: an antiproton beam momentum of 
1.5 GeV/ c and the detection of the lepton pair. The count 
rate in the near-threshold region of Ai is small but not 
negligible. The main difficulty is to identify the reaction 
among the hadronic background which is about six or- 
ders of magnitude higher. First elements of strategy were 
presented for this background rejection, based on parti- 
cle identification, detector hermeticity, and missing mass 
resolution. Although the subject would require a much 
more detailed study, we conclude that this process has a 
chance to be measurable in PANDA, given the very good 
design performances of the detector. 
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